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ABSTRACT 

When the gas of a magnetized filamentary cloud obeys a polytropic equation of state, 
gravitational collapse of the cloud is studied using a simplified model. We concen- 
trate on the radial distribution and restrict ourselves to the purely toroidal magnetic 
field. If the axial motions and poloidal magnetic fields are sufficiently weak, we could 
reasonably expect our solutions to be a good approximation. We show that while 
the filament experiences gravitational condensation and the density at the center in- 
creases, the toroidal flux-to-mass ratio remains constant. A series of spatial profiles 
of density, velocity and magnetic field for several values of the toroidal flux-to-mass 
ratio and the polytropic index, is obtained numerically and discussed. 
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1 INTRODUCTION 

Understanding the processes of transforming molecular clouds into stars is one of the main goals of the people who are working 
on the structure formation in interstellar medium (ISM) . Since various physical agents such as self-gravity, thermal processes, 
magnetic fields are playing significant roles in star formation, we are still far from a coherent and consistent picture in spite of 
great achievements during recent years. While the clouds in the standard model of star formation are considered to be initially 
in equilibrium or quasi-hydrostatic phase, some authors questions whether this important aspect is plausible according to the 
numerical simulations and the observations (e.g, Vazquez-Semadeni at al. 2004). Irrespective of which theory can correctly 
describe the initial state of the clouds, the next stage of evolution of the cloud is gravitational collapse. There are many studies 
for presenting a correct description of the gravitational collapse of the clouds, considering their geometrical shapes and the 
main physical factors. 

Studies of the gravitational collapse of gaseous clouds have been started by the pioneer works of Bodenheimer & Sweigart 
(1968), Larson (1969), and Penston (1969) who studied the isothermal collapse of spherical clouds using numerical integration 
of the equations and semi-analytical similarity solutions. It seems that self-similar fiows provide the basic physical insights to 
the gravitational collapse, and may indicate the way to more detailed investigations (Shu 1977; Hunter 1977, 1986; Whitworth 
& Summers 1985). 

Most of the previous works dedicated to study of the collapse of spherically symmetric clouds. But it is known that 
filamentary structures associated with clumps and cores are very common (e.g., Schneider & Elmegreen 1979; Houlahan & 
Scalo 1992; Harjunpaa et al. 1999). Also, filamentary structures are so prominent in numerical simulations of star formation 
that seems the formation and fragmentation of filaments is an important stage of star formation (see, e.g., Jappsen at al. 2004). 
So, self-similar solutions for the collapse of a filamentary cloud were investigated, and different sets have been found (Inutsuka 
& Miyama 1992; Kawachi & Hanawa 1998; Semelin, Sanchez, & de Vega 1999; Shadmehri & Ghanbari 2001; Hennebelle 2003; 
Tilley & Pudritz 2003 ). Recently, Hennebelle (2003) investigated self-similar collapse of an isothermal magnetized filamentary 
cloud which may undergo collapse in axial direction, in addition to radial direction. Then, Tilley & Pudritz (2003) extended 
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this analysis by focusing on only the purely radial motions and obtained interesting analytical solutions able to describe 
collapse of an isothermal magnetized filamentary cloud with purely toroidal magnetic field. 

While most of the authors assumed that the clouds arc isothermal and so their self-similar solutions describe collapse 
of isothermal filamentary clouds, Shadmehri & Ghanbari (2001) studied quasi-hydrostatic cooling flows in filamentary clouds 
using similarity method. Although an isothermal equation of state is a natural first approximation, there are some growing 
evidences that precise isothermality is not expected in molecular clouds. Scalo et al. (1998) studied the likely values of the 
exponent 7 that appears in polytropic equation of state of the form P = Kp'' . Kawachi & Hanawa (1998) investigated the 
gravitational collapse of a nonmagnetized filamentary cloud using zooming coordinates (Bouquet et al. 1985). They used a 
polytropic equation of state to indicate the effects of deviations from isothermality in the collapse. 

In this paper we derive self-similar solutions that can describe collapse of a polytropic filamentary cloud considering 
magnetic effects. While in the isothermal case the radial velocity is in proportion to the radial distance and one has to choose 
a prescription between magnetic field and the density (Tilley & Pudritz 2003; hereafter TP), we show that when the gas 
obeys a polytropic relation not only the radial velocity has not such simple nature but also the toroidal flux-to-mass ratio is 
constant during the collapse phase. The equations of the model are presented in the second section. We obtain and solve the 
set of self-similar equations in the third section. These solutions will be discussed in this section. 



2 GENERAL FORMULATION 

In order to study gravitational collapse of filamentary magnetized clouds, we start by writing the equations of ideal magneto- 
hydrodynamics in cylindrical coordinates (r, (p, z). We consider axisymmetric and long filament along the z axis. Thus, all the 
physical variables depend just on the radial distance r and time t. As for the magnetic field geometry, the toroidal component 
of the field is assumed to be dominant. The governing equations are the continuity, 

l + w 

the momentum equation, 

dv, dv, Idp 9* d , „ ^ 

15" +^■^"5" + + ^ -S-r-Bv): 

at or p or or ppr or 

the Poisson's equation. 

Id 9* 

-^ir^)=^^Gp, (2) 
r or or 

and the induction equation, 

f- + |:(..B.) = o, (3) 

We also assume a polytropic relation between the gas pressure and the density, 

p = Kp'', (4) 

with K and 7 = 1 + ^ are constants and 7 < 1 (Scalo et al. 1998). So, we relax the isothermal approximation adopted in 
the previous study by TP. Since we neglected the axial velocity Vz, the corresponding solutions could only be applied to the 
regions near the middle of a filament having finite length. 



3 SELF-SIMILAR SOLUTIONS 
3.1 similarity equations 

In the self-similar formulation, the various physical quantities are expressed as dimensionless functions of a similarity variable. 
The two-dimensional parameter of the problem are K and Newton's constant, G, from which we can construct a unique 
similarity varaible 

^ = K-iG'^r{to-ty-'\ (5) 

where t < to and the term, to, denotes an epoch at which the central density increases infinitely. Dimensionless density, 
velocity, gravitational potential, and toroidal component of magnetic field then be set up as 

p{r,t) = G-'{to-t)-^R{0, (6) 
p{r,t)=KG-'{to-t)-'-'P{0, (7) 
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Vr{r,t)^K-^G^{to-ty-''V{0, (8) 

B^{r,t) = ii^K^G-^ (to -tyh^ (0, (9) 

*(r,f) = KG'-^to - tf^'-'-'^SiO- (10) 

Thus, in terms of these similarity functions, the equation of state simply becomes P = W . The continuity equation, Euler 
equation, Poisson equation, and induction equation become 

2R+{2-^)i'^ + )^mV)=Q, (11) 
-,^.r /o -..dV ,,dV 1 dP 



76^ + (2-7)?^ + ^(V'6^)=0, (14) 

We can obtain solutions of TP for the collapse of an isothermal magnetized filament, simply by substituting 7 = 1 in the 
above equations. In this case, equation (11) is integrable and gives V = Thus, equation (14) is automatically satisfied and 
one has to choose a relationship between the magnetic field and the density in order to make further progress (TP). Clearly, 
in the polytropic case 7^1, the continuity equation (11) has not the simple nature of the isothermal collapse. However, from 
Equations (11) and (14) one can simply show that the toroidal component of magnetic field, 6<^, should be proportional to 
^_R. On the other hand, the toroidal fiux-to-mass ratio r,^ is defined as F,^ = B^/rp (Fiege & Pudritz 2000). Considering 
equations (5), (6) and (9), we can write F^ as 

r. = ||. (15) 

Since the above similarity equations show that h,p is in proportion to ^i?, the toroidal flux-to-mass ratio F,^ should be constant. 
We can consider F<^ as free parameter. 
If we write 

C^ = (2-7)^ + ^^, (16) 
then we can re-write the equations as 

dR _ RX , . 

D ' ^ ' 

dV {2 + V/OD + UX 

dT ~ D ' ^^^> 

where 

D = ■yR''-^ -U^+ T%i^R, (19) 

X = {2 + ^)U + {1-^)V - -^RU -2VliR. (20) 

? 1-7 

Equations (15), (17) and (18) describe the gravitational collapse of a polytropic magnetized filament. These equations are 
very similar to the equations obtained by Larson and Penston for the collapse of a spherical non-rotating and non-magnetized 
cloud (Larson 1969; Penston 1969) and extensively investigated by Hunter (1977), Shu (1977) and Whitworth & Summers 
(1985). If we set = 0, the equations reduce to the equations derived by Kawachi and Hanawa (1998) for the collapse of a 
polytropic unmagnetized filament. The aim of this paper is to solve these coupled differential equations (17) and (18) subject 
to suitable boundary condition, and interpret the solutions. When D vanishes, all the numerators in equations (17) and (18) 
must vanish at the same time; otherwise the equations become singular and will yield unphysical solutions. 

3.2 special analytical and asymptotic solutions 

Equations (17) and (18) have an analytical solution, 

y = 0, (21) 
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Figure 1. Profiles of the density, the radial velocity and the toroidal component of the magnetic field corresponding to 7 = 0.9 (top), 
0.8 (middle) and 0.7 (bottom). Each curve is marked by the toroidal flux-to-mass ratio T^. 



R = A£,~, (22) 

h^^V^AC^ , (23) 

where A is a constant. In dimensional units, this solution corresponds to time-independent singular configuration. 

We can find a family of limiting solutions for ^ ^ 0. In order to obtain these, one can expand -R(^) and V(£,) in Taylor 
series 

R{(.) ^Rq + + R2f + ■■■ (24) 

V(0 ^Vo + ViC + V2^^ + ... (25) 

By substituting the above series in equations (11)-(14), we obtain the following limiting solutions for ^ parameterized by 
Ro and F^: 

RiO -Ro- + + oi^), (26) 
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Table 1. Summary of the similarity solution for various 7 and F, 
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1.10 
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1.24 


731 


2.24 


1.46 


124 


9 


2.37 


0.78 
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1.25 
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2.41 
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2.36 
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10446 


2.35 


1.07 


1386 



nO = -C-^|^«o""^' + 0(C'), (27) 

= r^^iRo - ""^m'^ Rl-^e + o{t)]. (28) 

Asymptotic solutions (26)-(28) are very useful when performing numerical integrations to obtain similarity solutions starting 
from ^ — » 0. 



3.3 numerical solutions 

When similarity solutions cross the critical points, wo must take extra care as for their analyticity across the critical points. 
In order to find the behaviour of solutions around critical point, ^ = we expand R{£,) and V{^) in a Taylor series up to 
first order. Detailed calculations are in Appendix A. Given 7 and F^, we can find analytical expressions for ^« and V* as 
functions of 7?*. However, note that only those points are acceptable that ^, > and V* < 0. Although we don't study type 
of critical points in detail, when performing numerical integrations we note if the critical points are saddle or nodal (Jordan 
& Smith 1977). In general, the critical point is a saddle point when the two gradients of eigensolutions of V(^) are of the 
opposite signs. On the other hand, when the two gradients have the same sign, the critical point is a nodal point. 

We can numerically integrate equations (17) and (18) by a fourth order Runge-Kutta integrator. By using asymptotic 
solutions (26) and (27) for ^ — > 0, one can start integration of equations (17) and (18) at a very small ^ ( e.g., ^ = 10"*) 
and with arbitrary value for Ro for given and 7. Also, we can integrate the same equations (17) and (18) backward from 
the critical position, using asymptotic behaviour of solutions near the critical point for an arbitrary Thus, we look for 
matchings of R and V at some specific point (0 < < C*)- This can be found by studying loci in V — R plane at ^m. 

Table 1 summarizes the dependence of the similarity solutions on and 7. Listed axe the position of the critical point 
(^t), the density at this point (Rt.) and the density at the center {Ro). Figure 1 shows profiles of the density R, the radial 
velocity V and the toroidal component of the magnetic field 6^ for different choice of 7 and F^. Clearly the solutions except 
for the geometry and the magnetic field axe equivalent to the Larson-Penston solution. In almost all cases, there are an inner 
part with approximately flat density profile and an outer region with decreasing density. Since the toroidal fiux-to-mass ratio 
F(^ is constant during the collapse, each curve in Figure 1 is labeled by this ratio. For a fixed F^, the central density increases 
as 7 increases. The infall velocity increases and tends to an asymptotic value at large radii. Of course, such large velocities are 
not acceptable and there must be a mechanism (e.g., pressure turncation) to turncate the solutions at finite radius. However, 
the typical behaviour of the radial velocity of a collapsing polytropic filament is different from the isothermal collapse where 
the self-similar infall velocity behaves in proportion to the radial distance. Also, Figure 1 shows that in the inner part the 
toroidal component of the magnetic field increases, while in the outer part b^ is a decreasing function of ^. This behaviour 
is easily understood, if we note = F^^^Ti. For example, in the inner part, the density is roughly constant and so b,p oc 

An interesting feature of the solutions is that the cloud is more compressed by the toroidal pinching, if the toroidal 
magnetic field increases. Figure 1 shows that the size of the inner region strongly depends on F<^. However, the typical 
behaviours of the density and the magnetic field in the outer part are more or less independent of the toroidal field. We see 
that the central density increases, if F,^ increases so that in highly magnetized filament, the inner part with flat density profile 
disappears and joins to the outer region. If we compare the density profiles of 7 = 0.9 and 7 = 0.7 for F<^ = 10, we see there 
is still a small inner region for lower 7 = 0.7. It simply implies as 7 decreases, there needs higher level of magnetic intensity 
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in order to affect tiie ffat density profiie of the inner part. In other words, as the cloud tends to the isothermal regime, the 
inner part becomes more sensitive on the toroidal magnetic field. 

The infall velocity reduces as increases. While in nonmagnetic collapse, the radial velocity tends to very high value at 
large radii, we see that the velocity significantly reduces even at large radii in highly magnetized collapsing filament. However, 
the density profile in the outer region hardly depends on T^. We find this profile mostly depends on 7, so that it behaves 
in proportion to r^^'^^ , r^^''''^ and r^^"''' at large radii for 7 = 0.9, 0.8 and 0.7, respectively. This typical behaviour of the 
density is in good agreement with the asymptotic solution at large radii which is expressed as R{^) = Roo^~'^^^'^~''^ ■ Behaviour 
of the density profile at large radii does not change by increasing F,^ for a given 7. This behaviour is different from previous 
studios. Stodolkiewicz (1963) and Ostriker (1964) studied isothermal equilibrium structure of an isothermal filament, which 
the density falls off as r~* at large radii. Miyama et al. (1987) derived a set of self-similar solutions for an unmagnetized 
collapsing filament and their density structure is similar to unmagnetized equilibrium solutions, i.e. in proportion to r~'^ at 
large radii. TP showed that the density profile of an isothermal collapsing filament with purely toroidal magnetic field may 
change at large radii depending on the level of magnetization from (low magnetization) to r^^ (high magnetization). But 
we see as the cloud deviates from isothermality, these typical behaviours of the density profile at large radii change and the 
main parameter in shaping the profile is 7 not the level of magnetization, i.e. F^^. 

Considering profile of the density in outer part of the filament, we can obtain the toroidal component of the magnetic 
field. Since at large radii we have R oc r", the toroidal field becomes 6<p oc r""*"^, where v = —1.81, —1.69 and —1.55 for 7 = 0.9, 
0.8 and 0.7, respectively. We discussed that the density profile at large radii is fairly insensitive to F^, and it simply implies 
that the scaling of with the radial distance in the outer part is also insensitive to F<^. We can find behaviour of the ratio of 
the thermal to the magnetic pressures, f3. One can easily show that in our notation the ratio is /3 = {2/T'^){R'^~^ /(^^) . Thus, 
we see that in the inner part f3 oc r~^, irrespective of the value of 7 and F^. But in the outer region we have (3 oc r'^^'^'"^'^. 
These behaviours of /3 show that during gravitational collapse of a polytropic, magnetized filament this ratio is not constant. 
However, the toroidal flux-to-mass ratio is conserved during collapse phase. In isothermal regime, one can assume either 
/3 = constant or F,^ = constant and then study gravitational collapse (TP). 



4 CONCLUSION 

In this paper, we have derived solutions able to describe dynamics of a collapsing, polytropic, magnetized, self-gravitating 
filament. The magnetic field is assumed to be purely toroidal, although it is the only possible magnetic configuration which 
can be studied using the similarity method (sec below). But we expect our solutions to be a good approximation, if the 
poloidal fields are sufficiently weak. Moreover, purely toroidal field makes it easier to compare the obtained solutions with the 
isothermal collapse solutions as TP studied. Also, The self-similar solutions of our model can be used for future numerical or 
analytical studies because except for the geometry and magnetization they resemble to the Larson-Penston solution that has 
been found to be in good agreement with numerical analysis (e.g., Larson 1969; Hunter 1977; Foster & Chevalier 1993). 

One may ask is it possible to extend this analysis by considering both the poloidal and the toroidal components of 
the magnetic field. We note that fiux conservation imposes Bz oc p, which means the similarity solutions should support this 

11 T 

scaling. However, such similarity solutions should be in these forms, according to a dimensional analysis: Bz = /i^ 7^2 2 (£0 — 
£)~^6z(C) and p(r, t) = {to — 1)~^ R{£,) . Obviously, this scaling does not support the fiux conservation and so, it is necessary 
to put fez = or Bz = . On the other hand, Hennebelle (2002) explored a set of self-similar solutions for a magnetized 
filamentary cloud, in which the toroidal and the poloidal components of the magnetic field play significant role in the dynamical 
collapse of the cloud. However, his solutions describe an isothermal, magnetized filamentary cloud which undergo collapse in 
the axial direction, in addition to radial collapse. In his model, the slope of the axial velocity being two times the slope of the 
radial one at the origin. 

One of the major conclusions of our study is that the dynamics of a polytropic filament is different from the isothermal 
case. We showed that the toroidal component of the magnetic field help to confine the gas by hoop stress and this conclusion 
is independent of the exponent of the polytropic equation of state. Most measurements of molecular clouds have difficulty 
resolving the inner regions of filaments and so, it is very important to understand behaviour of the physical quantities at large 
radii. Our solutions showed that the typical behavior of the density of a magnetized polytropic filament mainly depends on 
polytrop index, 7, irrespective of the level of magnetization. However, the infall velocity in the outer regions strongly depends 
on the ratio of the flux- to-mass ratio . It implies to examine the profiles of molecular lines for infall velocity which may 
help us to understand the true nature of the filaments. While radial velocity of the isothermal model of TP has a very simple 
behaviour, our polytropic solutions and Hennebelle (2002) model clearly present different profiles of the radial velocity similar 
to Larson-Penston solutions of collapsing unmagnetized spherical clouds. 

The environment of L1512, a starless core, has been studied at high angular resolution by Falgarone, Pety & Philips 
(2001). The gas outside the dense core is structured in several filaments with a broad range of density and temperature, from 



© 0000 RAS, MNRAS 000, 000-000 



Gravitational collapse 7 



uhj = 2 X 10^ cm~^ for the coldest case (T = 20 K) down to nn^ = 180 cm~^ for the warmest (T = 250 K). It suggests that 
a polytropic equation of state is more appropriate for describing the thermal behaviour of filaments of this system. Falgarone 
et al. (2001) discussed that these filaments arc not held cither by the pressure of the H i layer or by the external pressure. 
So, they concluded that the toroidal component helps confine the filaments. However, their arguments were based on the 
analysis of Fiege & Pudritz (2000) , in which the equilibria of pressure-truncated isothermal and logatropic filaments held by 
self-gravity and helical magnetic fields have been studied. However, it is unlikely that the filaments in ISM arc being truly 
static structures, and many numerical simulations show dynamics structures in a star forming region (e.g., Vazquez-Semadeni 
et al. 2005; Jappsen et al. 2004). We think dynamics models, like what studied here, are more adequate for filamentary 
structures such as those in the environment of L1512 which show a wide range of density and temperature. 

Acknowledgements: I thank the referee, Anthony Whitworth, for a careful reading of the manuscript and comments that 
lead to improvement of the paper. 
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APPENDIX A: 

To obtain accurate transonic solution, it is useful to analyze the behavior of the fiow near the sonic point, ^, > 0. The values 
of = R{i*) and K = V(^*), where V* < 0, are completely determined by requiring both the denominator and numerator 
of equations (17) and (18) to vanish at 

7iiZ-' - + V%CR* = 0, (Al) 

(2 + ^)U, + (1 - 7)K - -^^RM* - 2r%^,R, = 0, (A2) 

where Ut = {2 — 7)^* + Vf In general, given 7?,, it is a simple matter to find ^« and K from equations (Al) and (A2). After 
mathematical manipulation, we obtain 



, _ , -{ca + /) ± Vica + ff + 4/i(c6 - d) ^,^, 
^* - ^ 2(c6-d) J 

where 

a = ^RZ-\h = VlR^.c = (1 - ,^R,f, 

d = [(1 - 7)(2 - 7) -h / = 2ad'/', /i = a\ 

© 0000 RAS, MNRAS 000, 000-000 



8 Mohsen Shadmehri 



Now, wc set 

7?(e)=i?*+ai(C-e*) + 0[(C-e.)'], 
V{0^V, + b-,{^-^,) + O[{^-^,)% 

in neighborhood of Substitution of this first-order expansion in equations (17) and (18) leads to the following equations 
for the slopes ai and bv 

aiDi - R^Xi = 0, (A4) 
{2 + bi + ^)Di+U.Xi=0, (A5) 
where 

J)[/. + (2 + ^ 



-^1 = (f- - 75-)^^* + (2 + ^)(2 - 7 + 6i) + (1 - 7)61- 



^^[oi?7. + i?»(2 - 7 + 61)] - 2r^(Ji. + ai^*), 

Di = 7(7 - l)Rr^ai - 2(2 - 7 + 6i)[7.+ 

r^C*(aie. +2i?.). 

Prom equation (A4), we have 

ai = -^(2 + fei + |l), (A6) 
and substituting in equation (A5) gives an algebraic equation for 61 as 

Abl + Bbi+C = 0, (A7) 
where 

A = -7(7 - 1):^ - 2U. - r^Ml, B = Bi+B2 + B3,C = U^Ci + C2) + (2 + ^)C3, 

B, = U4S-y+^ + ^ + 2r^^), B2 = (2 + |^)[^ - 7(7 - 1)^]' S3 = 24ii4* - 4(2 + ^)^ - 2(2 - 7)^/., 
a = (2 + |1)(2 - 7 + ^) - = 2r^(2 + ^)Ml _ 2[!^ + r^,]i?., Cs = Ba - 7(7 - 1)(2 + ^)^- 
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